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,__! Abstract 

The dynamical stability of a Julia-Zee solution in the AdS background in a four 
dimensional Einstein- Yang-Mills-Higgs theory is studied. We find that the model with 
a vanishing scalar field develops a non-zero value for the field at a certain critical 
temperature which corresponds to a topological dyon in the bulk and a topological 
phase transition at the boundary 



1 Introduction 

The gravity-gauge theory duality, AdS/CFT [I] [2] [3], has provided the hope of doing quan- 
titative analysis in low energy QCD. It is suspected that the phenomenon of confinement in 
QCD may be related to the condensation of topological configurations, specifically magnetic 
monopoles in QCD El [5] [6] . Extended topological configurations in field theories have 
other significant consequences. 

In this work we study the consequence of one such topological charged magnetic con- 
figuration, a dyon, in the bulk of an AdS/ CFT set up. The configuration being the dyon 
solution of an SU(2) gauge theory coupled to a triplet of scalar fields, of 'tHooft-Polyakov- 
Julia-Zee type, in the asymptotically AdS Reissner-Nordstrom blackhole background. We 
will argue that this will emerge in the low temperature regime from a configuration with no 
scalar field. The boundary theory will make a corresponding phase transition. 

Magnetic monopoles have a long history. The gauge theory in 4 dimensions, the electro- 
dynamics, does not allow magnetic charges at the right hand side of the Maxwell equation for 
the dual field intensity vanishes identically, unless a singularity, the Dirac string, is inserted. 
A similar situation in non-abelian gauge theories persists. 

Nearly thirty five years ago Polyakov and 'tHooft observed that introduction of an ad- 
joint scalar permits the appearance of a magnetic monopole solution in the SU(2) gauge 
theory, upon spontaneous breaking of the gauge symmetry to a U(l) subgroup determined 
by the direction of the scalar field in isospin space. The magnetic charge being a topological 
quantity, is classically quantized [7j [8J. 

This "Hedge-hog" solution led to an avalanche of work on generalization to other gauge 
groups with serious problems for cosmology, whose resolution was achieved in the theory of 
inflation. 

Shortly after the discovery of the solutions with magnetic charge, dyons with electric 
charge also were found by Julia and Zee J9J. A particular limit in the solutions discovered 
by Prasad and Sommerfield [10] , saturating a bound observed by Bogomol'nyi fllLthe BPS 
solution has permeated the literature of supergravity and string theory. 

Magnetic monopole solutions in the presence of gravity have also attracted a great deal 
of attention. Many solutions in asymptotically flat, deSitter, or Anti-deSitter spaces with 



and without blackhole singularity have been discovered 12 . Remarkably only the ones with 



the Anti-deSitter asymptotics are stable 13 14 . This is a striking result in view of the 
AdS/CFT duality. 

In particular for SU(2) gauge theory, asymptotically AdS blackhole solutions with and 
without scalar field have been discovered (15), and more recently by Lugo, Moreno, and 
Schaposnik, [l6| . 

We will use a particularly simple solution found quite some time ago by Kasuya and 



Kamata 17 , in the present work. We will find that their solution in the absence of a scalar 
field is unstable while the one with non-vanishing scalar triplet is. 

The paper is organized as follows. In the next section we briefly review the Julia-Zee 
dyon in flat space. In section 3 we introduce it's extension to gravitational versions. In 
section 4 we review the general procedure for studying dynamical stability. In section 5 we 
apply this procedure to our case and derive the main analytic equations. In section 6 we 
present the results of our numerical calculations. Section 7 is devoted to conclusion. 



When this work was completed, a paper by Lugo, Moreno, and Schaposnik appeared [18], 
with similar considerations, related to their previous solution and in the context of the 
noncompact boundary of R 2 in place of S 2 . 

2 Julia Zee dyon in flat space 



In 1974 'tHooft |7| and Polyakov |8| independently introduced a new type of magnetic 
monopole in a flat Minkowski space. Their monopole is free of any singularities. The 
action of the theory is Yang-Mills-Higgs with a special nontrivial ansatz for the gauge and 
scalar fields. The action is: 



S = J ' d 4 x i—F^F^ -i/VWT - V((f> a n 



(2.1) 



The scalar potential is assumed to have a minimum at <fr a <j) a = (0) 2 = const ^ 0. The 
covariant derivative and the field strength are defined as usual, 

D^ a = d,r + ee abc A b ^ c ; F^ = d,A a u - d„A; + ee^A^. (2.2) 

The 'tHooft-Polyakov ansatz for a solution for this system is, 

$ = I ^lf ; A % = i— ^a 4 (i = 1,2,3) ; < = e mk x k ; A t = (2.3) 

er er 

Regularity and finiteness of energy requires the boundary conditions, 

atr-^0 : F -» ; if ->• 1 

at r — )• oo : H — y const x r ; K — > 

The U(l) field strength, invariant under the su(2) algebra, is defined as, 

F^ = F^4 + l[D^ x Dj] ; j> = -$= ■ (2.4) 



When r — > oo, the magnetic field scales as 1/r 2 ; so the configuration has a magnetic charge: 

lim Bi ->• — ? =>- <p Bi ds i = - = g (2.5) 



r— >oo e r 2 I e 



Later Arafune, Freund, and Goebel 19 showed that this magnetic charge is a topological 
object and its value is discrete, 

g = i Bids 1 = -. (2.6) 

The hamiltonian of the system can be written as, 

H = Jd 3 x \^{F^-e l3k D k <p a ) 2 + V(r<P a )]+^g(<P), (2.7) 



where (0) is the boundary value of which is a constant by assumption. It is then easily 



found 11 that there is a bound on the energy, BPS bound, 

E > 4vr#(0). (2.8) 

The bound is saturated when 

Ft 3 = e ijk D k <f) a . (2.9) 

This magnetic monopole solution was later extended to a dyon by Julia and Zee [9]. 
Their solution was in the form: 



<j)= r, 

er 

K = ^^k, (2-10) 



A J ( r ) ~ 

A t = ——t-r. 
er 

The boundary conditions for H and K are as before, and the boundary conditions for 
J(r) are: 

at r ->■ : J —> 

at r — y oo : J — >• const x r 

With this ansatz we see that the configuration has an electric charge too. 

In 1975 Prasad and Sommerfield flO] found an exact analytic solution for this dyon in 
the limit V — y 0. their solution has the form: 

K Cr 



sinh{Cr) ' 

J = sinh{i) [Cr coth(Cr) - l] , ( 2 ' U ) 

H = cosh(j) [Cr cothiCr) - l], 

where C and 7 are arbitrary constants. 

3 Julia- Zee dyon coupled to gravity 

The dyon of the previous section lives in flat space but what happens if we include the 
gravity? During 80's it was shown that there exist gravitational Julia-Zee dyons; and certain 



exact solutions were obtained 20 . One of the simplest solutions in this class is the Kasuya 
and Kamata solution (l7|. 



When gravity is included the action is: 

s = Jv^gd'x [^q{r + jp)- \f; u f^ - l -D^ a D^ a - x v(rr)} (3.1) 



And as before the scalar potential is assumed to have a minimum at some nonzero constant 
(4>) 2 . The equations of motion are: 

1 6 

-^V — 9^ + 72) 9nv = ^G T^, 

T»„ = [^ - \F; x F a ^ x - ^D x rD x r - A\/(0 2 )]^ + f; x f: x + D^ a D„r, 

d,(^-gD^ a ) + V^g{ee abc 4Dy c - A ^) = 0, 

d»(V^9F a ^) - eV^g e abc (F b ^A c ^ + (<9> 6 ) C ) + e 2 v^((A 6 > fe )0 a - (<j> b <i> b )A av ) = 0. 

(3.2) 

The same ansatz for scalar and vector fields is assumed as in the flat space case. The 
most general form for a spherically symmetric metric in 4 dimensions is: 

9ilv = diag(-e x{r \ e r(r) , r 2 , r 2 stn(6) 2 ) (3.3) 

If the scalar potential is taken to be in the form of a mexican hat, V = |(0 a a ) 2 — \{4> a (j) a )- l 
then the equtions of motion for H, K and J become: 

r-Ux + ¥)'£)' -e Y ^ = o, 

2 r r z 

K" + -(X - Y)'K' - e Y - 2 {K 2 + H 2 - e~ x J 2 - l) = 0, (3.4) 

H" + r -(X - y)'(f )' - e Y y 2 (2K 2 + ±(H 2 - cV)) = 0, 

where c 2 = e 2 v 2 and v is the minimum of the scalar potential. 
The equations of motion for the metric yield, 

-Y 



e~ Y , „, ., 1 8ttG 



(rX' + l) + 4 = ^T;, (3.5) 



/yiZ /yiZ zD -i 



■^ [a* + i(A-<) 2 - irv" + V - Y)>] = ^T" = ^7J, 



r e^ e^ 



where Tf,T^,Tg and T£ are the components of the energy momentum tensor, 

-Y 



T 



iK , ) 2 + K -^-T^ + —^[(-)T+ 



7-2^2 



x J 2 K- e 

1 ( — 

r 4 2 [K r 



(K 2 - l) 2 e-( x+y ) r/ J v1 
2T 4 + 2 
" ff vi2 , tf 2 ir 2 , , 



+ e 2 Al/(0) 



rpr 



-y 



■{K'f + 



{K 2 



1) 



2^2 



-x 



J 2 K 



-y 



H, 



-K 



/12 



2r 4 
# 2 if 2 



+ 



.->■■» 'j 



-(x+y) t 



+ e 2 \V{(t>) 



(3.6) 



7^6 



rpip 
<P 



(K 2 -l) 2 e~( x+Y \J v 



-Y 



2r 4 



[(■ 



+ 



,H 



— [(-)'Y + e 2 XV{<P) 



Kasuya and Kamata found an exact solution to these equations with the above ansatz 17 

,2 



where, 



H — cr ; K — ; J = \xr — p \ e x = e Y 

2 _ 47rG w (l + p 2 ; 



2M o 2 r~ 

1 +^+ V7 



r 



<f 



(3.7) 



(3.J 



This expression is still a solution if we set H = 0. The background metric in both cases 
are AdS-RN; but, there is an important point : The AdS radius differs for the two, because 
the contribution of the scalar potential to the cosmological constant is different for the two 
cases. Note also that because the gauge field should have a finite norm, A\ should be zero 
at the horizon of the blackhole, implying a relation between \i and p 



I' 



_P_ 



(3.9) 



where r# is the outer horizon of the blackhole. The temperature of the blackhole after 
eliminating the charge of the blackhole from definition of the horizon (e x = 0) is: 



1 , M 2r 2 H , 

T = (l 1 -). 

2nrH Tn L 2 



(3.10) 



Concerning the Bogomol'nyi equation, when the gravity is included, the simplest guess 



is: 



D: 



-g £ijk 



pajk 



(3.11) 



But this does not work; in fact this relation is not compatible with the equations of motion 
derived from the action (3.1). It has been shown that the action of the theory can be changed 
in such a way that its equations of motion are compatible with the relation (3.11) |21|. On 
the other hand in certain circumstances a Bogomol'nyi like equation can be written which 



is compatible with the action (3.1). For example one may consider a generalization of the 
Bogomol'nyi equation (2.9) of the form, 



A0 a + di(u) r = ^V^9 e ijk F a ^\ (3.12) 

where u is an additional function. Consistency of this equation with the equations of motion 
leads to: 

u = log(y/\gm\ ) ; Aw = 0, (3.13) 

where A is the usual covariant Laplacian; this in fact is a constraint on the metric. For more 
details see |22|. 



4 General aspects of dynamical stability 

In this section we review the general aspects of dynamical stability analysis and holographic 
phase transition. 

In the presence of a background solution for a system, checking dynamical stability is to 
study the time evolution of the solution. Consider a system with a set of fields Xi(x) and 
the equations of motion Ei[Xj(x)] = 0, and a set of exact solutions X\ (x). Varying the 
solution with an infinitesimal time dependent variation, the linearized equation will be, 

L i [xf i) ,6X j (x),u]=0, (4.1) 

with SXi(x,t) = 8Xi(x) e~ iwt . 

If these equations admit a solution with an u which has positive imaginary part, then 
the variations blow up with time and the background solution will be unstable. 

This simple argument is the core of the concept of dynamical instability, but there are 
a number of technical issues concerning the boundary conditions on the fields which are 
important: 

The variations should be regular at any point of the space, as we assumed these variations 
to be infinitesimal. 

Also, at the horizon because of the general properties of the non-extremal blackholes, the 
radial dependence of the variations near the horizon is generally of the form [23] , 

6Xi(r) = {r-r H ) ±i ^r(a + a l {r-r H ) + a 2 {r-r H ) 2 + ...), (4.2) 

where T is the temperature of the blackhole, and the ± signs indicate the ingoing and 
outcoming modes. Then as a classical blackhole devours everything, the solution near the 
horizon should be ingoing only. 

Moreover, at the boundary at infinity, because the equations are of order two, the solution 
asymptotically has the form, 

5X ~ d r A + + c 2 r A " . (4.3) 

There are then two possibilities: 



1) One of the modes diverges (non-normalizable) ; the diverging mode must be excluded 
as the variation should remain small; while the other one vanishes at the boundary (normal- 
izable) . 

2) Both modes vanish. Then the one that is coupled to the appropriate boundary operator 
should be chosen. In the context of AdS/CFT such normalizable modes are interpreted as 
the v.e.v of the corresponding operator at the field theory (CFT) side. 

So our boundary conditions are: 

near horizon — > ingoing mode, 

(4.4) 
at the boundary — > C\ r A+ + c 2 r A ~. 

The coefficients c\ and c 2 are functions of parameters in the theory, e.g. temperature, radius 
of the blackhole, and chemical potential. Then choosing one of the boundary modes, e.g. 
c\ = 0, leads to the allowed frequencies: 

d(T,u,r H ,fj., ...) = ->■ u(T,r H ,fx,...). (4.5) 

Then if the imaginary part of the u changes sign in the vicinity of a hypersurface (wall 
of marginal stability) in the ( T,r H , /i, ...) space, then there will be a stability /instability 
transition for the configuration, and the system goes through a phase transition. 

The modes on the wall of marginal stability with u equal to zero are called marginally 
stable modes [24l. 



5 Stability analysis of gravitational Julia-Zee like so- 
lution 

We consider the quartic potential: 

\v(rr) = \(<P a <P a ? - ^(W- (s.i) 

As mentioned in section 3, there are two exact solutions to the equations of motion for 
Julia-Zee ansatz in the form of (2.10): H = with AdS radius L; and H = r yffi with 
AdS radius L : 

L-^ + 8 _^E f5 2) 

L2 " L2 + 3Ai> ' [b ' 



Perturbing the matter fields, 



er er e 

-^- + ee— '^, (5.3) 

er er e 

e e e 



where e is an infinitesimal parameter, and putting the new fields in the equations of motion 
for the solution H = 0, to first order in e, we get: 

o 
equation: (r 2 e x f ')'+(— + e~ x iu 2 ) r 2 f = 0, 

A t equation : 2P' + r P" = 0, 

A r equation : P' — 0, 

.4, e ,„ «on: { «'™ + e-[^ + (,-pAfW]Q = 0, (5 ' 4) 

A, equation: { V = + W + «-»[«*/.* + 0.-p/r)»+ ^9 = 0, 

The solutions to the gauge field perturbations are simply P = const , Q = 0; but, the 
equation for the variation of <p is nontrivial. We can consider different perturbations to A v 
and Ag, but the A t equation implies that they should be equal. 

The first step is to determine the boundary behavior of the function /. The near horizon 
equation for / is: 

f'+ f - + C- + %)f = 0, (5.5) 

u u u^' 

where u is (r — th) , and, 

' ^ = 7^, (5-6) 



' 2ttTL 2 ' 4ttT' 

with the solution, 

•W^V^), (5-7) 

which has the small u expansion, 

Z.±ivCj2 £.l±iv£p £,2±i\filr* 

u±lV ^ , _Jl V u _n ^ 3 y 

y r[l±2iV&] r[2±2iv^] 2r[3±2iV^] 

This corresponds to the form (4.2). 

The boundary equation for / at r — > oo is : 



nr» \ rr> Zi ryQ 



Changing the variable to z — 1/r, 



r-^r+(| + iv)/=o, (5.io) 

which has the solution, 

f = ze ±lL ^ z = -. (5.11) 

r 

For large r this has the form, 

/~7 + ^ + - (5- 12 ) 



As described in the last section a choice has to be made; at the horizon the ingoing mode 
is to be chosen and at the boundary at infinity either f\ or f 2 should be set equal to zero. 
Each leads to a different value for the critical temperature. 

We can repeat exactly the same procedure for the second solution with H — r yffi- 
At the linear order the equations are as before except for the <fi equation, 

(r 2 e x f>)'+(-^ + e- x ^)rZf = 0. (5.13) 

The forms of this equation near the horizon and the boundary at infinity are: 



11 11 11 



(5.14) 



where 

1 u L 1 . 

V = — i w = ~r^ i ot — — — — . , < a < 1, (5.15) 

' irTL 2 ' 4vrT ' L ^/TTj^ ~ ~ l ' 

7 is (87rGAr)/(3A), u = r — r^, and " ' " in the first equation denotes derivation with respect 
to u and in the second equation with respect to r. 
The solutions to the near horizon equation are: 



W>*( 2 V«$> ( 5 - 16 ) 

which are ingoing and outcoming modes. The solutions to the equation at the boundary at 
infinity are: 

1 ,La 2 u. , . 

^3/2 J ± | V9+16 Q 2 ( —^7~ ) , (5.17) 

with the asymptotic expansion of the form, 

/ = fir A+ + / 2 r A -, 

. -3 ± ^9 + 16a 2 ( 5 - 18 ) 

± 2 

The first term is divergent and is not acceptable, thus we set j\ = 0. 

6 Numerics 

We follow the procedure outlined above for finding the phase transition temperature by 
setting u equal to zero, i.e., finding the marginally stable mode. 



6.1 The case = 0: 

6.1.1 marginally stable modes: 

The equation when u> is set to zero is: 

(r 2 e X /')' + ^ r 2 / = 0. (6.1) 

The forms of this equation at the horizon and at the boundary at infinity are: 

r+- + -/=o, 

4 M 2 M (6-2) 

/" + -/' + - 2 f = 0, 

where u = r — r #; and " ' " denotes derivation with respect to u in the first equation and to 
r in the second one. 

The solutions to the first equation are: 



J (2^fju); Y (2y/rju). (6.3) 

Now Yo{2\/rpu) is divergent at the horizon so it is ruled out; thus we choose Jo(2-07 u), and 
its near horizon expansion is: 

1 -fju + .... (6.4) 

We can fix / at the horizon to be equal to 1. 
The solution to the second equation is, 

f = - + - 2 - (6.5) 

Eliminating the charge of the blackhole in favor of M, L and r#, then f\ and f 2 will be 
functions of (M,L,rn)- Imposing the desired boundary condition, f\ = or f 2 = 0, we 
will find a relation between (M, L, rn)- Simultaneously we have the positivity condition 
of the temperature , T(M, L, rn) > 0. For convenience we fix th to a certain value, e.g. 
th = 10, compute f\ or f 2 numerically and plot it as a function of M and L (by considering 
the positivity condition of blackhole temperature, —ML 2 + 2r# + r#L 2 > 0). Then we will 
have a surface, in the space of (M,L,f\) or (M,L,f 2 ) . Where this surface intersects with 
the plane f± — or f 2 — 0, is the critical curve of (M c , L c ). If there is no intersection then 
there is no marginally stable mode for the considered value of r# with the desired boundary 
condition at infinity. Our calculations show that for sufficiently large radius of the blackhole 
such an intersection exists. 

In figure (1) we have shown focused plots of f\ for two values of r#. For r# = 10 there 
is no intersection, but, for r# = 100 we have an intersection. 
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Figure 1: The green surfaces are the focused plots of f\ as a function of M and L and the 
blue planes are the planes for which /i vanishes ( plot (a) for r# = 10 and plot (b) for 
r H = 100). The plots are cut by T > condition. 



It can be seen that /i and f 2 can vanish several times before the temperature becomes 



zero. 



6.1.2 Quasi-normal modes: 

In this section we introduce our numerical results for quasi-normal modes. Our aim is to 
solve the first equation of (5.4) with the desired boundary conditions and find the complex 
frequencies. Our numerical calculations show that at high temperatures the = back- 
ground is stable, the imaginary part of the quasi-normal frequencies are negative; and at 
sufficiently low temperatures this background is unstable, the imaginary part of the quasi- 
normal frequencies are positive. 

It is obvious that we can not sweep the entire complex frequency plane by numerical 
calculations but it seems that at sufficiently high temperatures there is only one stable mode 
for every given temperature (or at least the frequencies in complex frequency plane are so 
widely distributed that we could only find one of them); but, at lower temperatures there are 
a number of unstable modes for every value of temperature. The quasi-normal frequencies 
that we found are listed in the tables below for both /i = and f 2 = 0. 

Note that, specially for /i = 0, in a range of temperatures above the onset of instability, 
the quasi-normal frequencies are purely imaginary or their real part is very small; this may 
be the sign of the existence of a mass gap right above the critical temperature T c . 

Note also that for f\ = 0, although there is no marginally stable mode at r# = 10, as 
mentioned in previous subsection, the quasi-normal frequencies exhibit a phase transition. 
This means that for the boundary conditions under consideration, when the temperature is 
changed, the frequencies do not go through the origin of the complex frequency plane. 
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Table 1: quasi- normal modes for /i = 



r H ==10 , L = l 


M 


4vr T 


U) 


2000 


0.212 


0.764969 + 1.66858 i 


2000 


0.212 


4.93743 + 3.35654 i 


1970 


0.812 


10.133 + 36.3126 i 


1970 


0.812 


0.669683 + 9.05345 i 


1950 


1.212 


9.9903 + 34.2773 i 


1900 


2.211 


0.0148391 + 35.9313 i 


1850 


3.211 


12.5418 + 52.0683 i 


1800 


4.210 


4.98854 x 10™ 15 - 1.30046 i 


1700 


6.210 


5.26857 x 10~ 14 - 2.09467 i 


1600 


8.209 


2.36149 x lO" 14 - 3.01488 i 


1500 


10.21 


1.26155 x 10~ 9 - 4.13518 i 


1200 


16.21 


2.16507 - 7.78775 i 


1200 


16.21 


8.66376 - 7.8904 i 


1000 


20.20 


3.7928 - 9.49205 i 



Table 2: quasi-normal modes for f% — 



r H = 10 , L = 1 


M 


4ttT 


UJ 


2000 


0.212 


1.06157 + 1.95297 i 


1990 


0.412 


1.09226 + 3.81241 i 


1980 


0.612 


1.0839 + 5.7522 i 


1950 


1.212 


3.31773 + 19.6673 i 


1900 


2.211 


1.77525 x 10~ 6 + 0.320354 i 


1500 


10.21 


0.0000327886 - 0.968159 i 


1400 


12.207 


0.0000337517 - 1.45497 i 


1000 


20.20 


1.55273 - 5.85023 i 


800 


24.20 


4.1098 - 5.3992 i 



6.2 The case 



2 . 
XL 2 ' 



We have also examined marginally stable modes for The Kasuya-Kamata solution (3.7). In 
this case the equations for / when ui is set to zero, are, 



full equation : 



(r 2 e X f)'-i- 2 r*f = 



L 2 



near horizon equation : f" -\ / = 0, 

u u 

4 ; Aa 2 
boundary at oo equation : f -\ — /' — / = 0. 



(6.6) 
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The solutions to the near horizon equation are: 



I (2^fju); K (2^fju 



(6.7) 



The second solution is divergent at the horizon so we choose the first one, which has the 
expansion : 

1 + fju + .... (6.8) 

The solution to the boundary equation at infinity is: 



/ - 



JiA + / 2 r A ". 
-3 ± V9 + 16a 2 



(6.9) 



The first term is again not acceptable so the second term has to be chosen. 

Finally we should find the appropriate set of parameters ( temperature ) for which the full 
equation admits a solution with the desired boundary conditions; and follow the procedure 
of the previous section, fixing the horizon condition and searching for parameters for which 
/ vanishes at the boundary at infinity. 

In figure (2) we show the plot of / at large r. We find that it does not vanish anywhere; 
therefore the equation does not admit a solution with the desired boundary condition at 
any temperature. Thus when (0) 2 = const ^ 0, the dyon solution, has no marginally stable 
mode, indicating stability, in accordance to its topological nature. 




Figure 2: The green surface is the plot of / at very large r as a function of M and L; and 
the blue plane is where f vanishes. The plot is cut by T > and q 2 > conditions. 



7 Conclusion 

In this work we studied stability of two simple dyon solutions of SU(2) gauge theory in an 
asymptotically AdS blackhole background and found that the solution without the scalar 
field is unstable. In the context of holography this is a phase transition for boundary theory. 
The solution with the non- vanishing scalar, describing a topologically nontrivial dyon, turned 
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out to be stable as expected. The SU(2) color symmetry and the SU(2) space rotational 
symmetry of the model is broken to a diagonal SU(2) symmetry in both phases. However 
the phase transition has a topological character. 

It is tempting to relate this holographic set up to a gravitational dual of a strong in- 
teraction physics model in 2+1 dimensions; the scalar iso-vector field of the model as pions 
which are believed to be Goldstone bosons of the spontaneously broken flavor symmetry. 
The color-spin locking in the dyon configuration is reminiscent of the color-flavor locking in 
QCD which is related to color superconductivity in neutron stars. And of course the intrigu- 
ing possibility of associating the topological magnetic configuration with the phenomenon of 
confinement in QCD. 

It would require further detailed study of the model to make any firm statements on 
these issues and we hope to get back to them. 
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